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A generalized Hamiltonian definition of ghost surfaces (surfaces defined by an action-gradient flow) is given 
and specialized to the usual Lagrangian definition. Numerical calculations show uncorrected quadratic-flux- 
minimizing (QFMin) and Lagrangian ghost surfaces give very similar results for a chaotic magnetic field weakly 
perturbed from an integrable case in action-angle coordinates, described by L = Lq + eL\, where Lq{6) (with 
6 denoting d6/d() is an integrable field-line Lagrangian and e is a perturbation parameter. This is explained 
using a perturbative construction of the auxiliary poloidal angle © that corrects QFMin surfaces so they are also 
ghost surfaces. The difference between the corrected and uncorrected surfaces is 0(e 2 ), explaining the observed 
smallness of this difference. An alternative definition of ghost surfaces is also introduced, based on an action- 
gradient flow in 0, which appears to have superior properties when unified with QFMin surfaces. 
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1. Introduction 

Recent calculations [ 1 1 of heat diffusion along chaotic 
field lines show that the isotherms correspond very closely 
with the "approximate" magnetic surfaces, associated with 
magnetic island chains, known as ghost surfaces [2 |. These 
surfaces include the "X-point" and "O-point" closed field 
lines of their associated islands. (By "O-point" field line 
we mean either the elliptically stable field line at the center 
of an island or its hyperbolically unstable continuation if it 
has undergone a period-doubling bifurcation.) Closed field 
lines extremize the magnetic action, j> A-dl, the hyperbolic 
X-point field lines in the chaotic separatrices being minima 
and the O-point field lines being minimax or saddle points 
of the action. Ghost surfaces are constructed by interpolat- 
ing smoothly between these two closed-field-line classes 
by evolving the O-point field lines into the X-point field 
lines along paths of steepest descent of action, thus gen- 
erating a family of "pseudo-orbits," i.e. paths that come 
close to extremizing action. 

An alternative approach to defining approximate mag- 
netic surfaces passing through magnetic islands, is to 
use the quadratic-flux-minimizing (QFMin) surfaces intro- 
duced by Dewar, Hudson and Price |2]|3]. Ghost surfaces 
have nice mathematical properties but are difficult to con- 
struct and have no obvious physical interpretation. QFMin 
surfaces, on the other hand, have the computational attrac- 
tion of being easy to construct, and the physical attraction 
of being defined in terms of a measure of the magnetic flux 
transport through the surface, but have been found to ex- 
hibit undesirable distortions in some circumstances. 

However the definition of QFMin surfaces is not 
unique, as it depends on the choice of an auxiliary poloidal 
angle 0; which raises the question: Can be chosen so 
that QFMin surfaces coincide with ghost surfaces? This 
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question has recently [4] been answered in the affirmative, 
provided we modify the definition of ghost surface slightly, 
there being some freedom also in the definition of action- 
gradient flow. This raises the further question, accepting 
that ghost surfaces need to be redefined to achieve unifica- 
tion with QFMin surfaces, is the definition used in Ref. |4] 
optimal, or is there a still better class of unified ghost and 
QFMin surfaces? 

In Sees. [2}]4]we review the QFMin and ghost surface 
concepts, and in Sec. [^summarize the unification proposed 
in Ref. 4 In Sec. [6] we summarize its perturbative imple- 
mentation in a test problem. 

In Sec.[7]we propose a new definition of ghost surfaces 
based on an action-gradient flow in the auxiliary angle 0, 
rather than the flow in the unperturbed canonical coordi- 
nate 9 used in Ref. [4] The new definition is not only sim- 
pler than that in Ref. 4 in that it does not require a variable 
rate of flow along pseudo field lines, it also leads naturally 
to straight pseudo-field-line orbits in the 0-£ plane. This is 
illustrated perturbatively for the same case used in Sec. [6] 




Fig. 1 A sketch of the general curvilinear toroidal coordinate 
system described in the text. 



2. Coordinates and Fluxes 

We use a general curvilinear coordinate system t/f, 6, 
as shown in Fig.[T] where 6 and g are poloidal and toroidal 
angles, respectively, and if/ is a label for the coordinate sur- 
faces that has the dimension of magnetic flux and is such 
that B = V^XV0 + V^xVy Q 

The standard linear flux JdS n-B through surface F: 
if/ = i//r(6, f)» where dS is an element of surface area and 
n is the unit normal, is independent of choice of coordi- 
nates and vanishes identically. Thus it is independent of 
the choice of F also, whether it be a magnetic surface (in- 
variant torus of the field-line flow) or otherwise. 

k= 1 
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Fig. 2 Poincare plot for the Standard Map with the p, q = 
1,3, 1,2 and 2,3 island chains, with associated action- 
minimizing (blue dots) and -minimax (red dots) periodic 
orbits. Also shown are almost-invariant curves interpo- 
lating between them. 



To measure the amount by which F departs from being 
a magnetic surface, we are instead led to define the positive 
definite quadratic flux |3] 
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where is an auxiliary poloidal angle, not necessarily 
the coordinate angle 0. Setting 5<f2 = for arbitrary 
variations i// r we find the Euler-Lagrange equation for 
QFMin pseudo-orbits pip) is B v -Vv = 0, where v = 
n-B/n-V0xV£ and the "pseudo magnetic field" B v is de- 
fined by 



B,iB - vV0xV^. 



(2) 



Thus v parametrizes a continuous family of pseudo field 
lines, which can be thought of as orbits of a dynamical 
system. We define a p, ^-periodic pseudo-orbit (p and q 
being integers) as a path that is a closed loop, with 6 in- 
creasing by 2np when ( increases by 2nq, so the average 



'Note that^ = x(i//, 8, cannot in general be made a surface function 
in three-dimensional systems due to the generic nonintegrability of the 
magnetic field-line flow. 



rate of increase of 8 along the path is the rational rotational 
transform t Pt9 - p/q. 

For a given v there are generally two distinct p, q- 
periodic orbits. In particular, when v = we find the true 
periodic orbits associated with the X- and O-points of a 
p, q island chain. 

The family of p, ^-periodic QFMin pseudo-orbits 
parametrized by v defines an "almost invariant" surface 
with rotational transform i p „, which passes through the 
p, q island chain as depicted in Fig. [2] 

3. Action and ghost surfaces 

Using the vector potential representation B = VxA, 
the magnetic action of a closed p, g-periodic path C is 



S[C] 



= A-dls A-tdC, 
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where r = dr/d£ along the path. Hamilton's Principle, 
6S = 0, gives the Euler-Lagrange equation 5S/dr = 
rxB = 0, which implies r is parallel to B as required. 

Closed paths making >S stationary are isolated, but the 
Poincare-Birkhoff theorem implies there are at least two 
associated with each island chain, as illustrated in Fig. [2] 
with the X-point orbits being minima of S, while the O- 
point orbits are minimax (saddle) points of S. The ghost- 
curve strategy for "joining the dots" is to generate a family 
of pseudo-orbits, labelled by a continuous parameter r, by 
flowing down the action gradient from minimax orbits to 
minimizing orbits, 
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where P e host is a symmetric nonnegative dyadic whose 
choice is discussed below. The union of this set of pseudo- 
orbits defines a ghost surface (a ghost curve being the 
Poincare section of a ghost surface). 




Fig. 3 A comparison of uncorrected (© = Q) QFMin curves 
(thick lines) and ghost curves (thin lines) for the more 
strongly chaotic case described in Ref. 4 Some cases 
where QFMin curves violate the graph property are seen. 



Ghost curves and uncorrected QFMin curves (i.e. 
curves for which the auxiliary angle has been taken 
to be the same as the coordinate angle 8) are compared 
in Fig. [3] For the lower-order p, q cases the two defini- 
tions of almost-invariant curves are almost indistinguish- 
able. However, higher-order uncorrected QFMin curves 
can cease to be graphs over 8 and do not agree with the 
isotherms found by Hudson and Breslau Thus we are 
led to seek a way to choose so that the two approaches 
may be reconciled, keeping the physical appeal and ease 
of implementation of the QFMin approach while retaining 
the superior mathematical properties of ghost surfaces. 

4. Lagrangian approach 

We first map this problem onto one of Lagrangian dy- 
namics in order to make better contact with the mathemat- 
ical literature, e.g. Ref.[5] 

The vector potential A = ifrVO - x(ip, 6, O^i gives 
our assumed form for B and leads to the action integral 
S[C] = f(iffd8 - xd£), the first variation being 

Xlnq 
[(6 - d t x) <ty - (<A + d eX ) 56\ d{ . (5) 

Setting 6S — for all diff and 68 leads to Hamiltonian 
equations of motion with 6 as generalized coordinate, if/ 
as canonical momentum, x as Hamiltonian, and f as time: 
6 = dyX and ip = -d x- 

The transition to the Lagrangian approach is made in 
the usual way, by solving 6 = d&x to give iff as a function of 
6 and writing the action as S — fLd£, where L(6,8,£) = 
tf/(6, 6 ,£)8 - xift, 8, <f ). The first variation now becomes 
5S - ^" c ' '56 (6S 1 '66) d(, the functional derivative (action 
gradient) being 

5S _8L d 8L 

~68 ~ 36 ~ d^dl ' 
Specialization of the general ghost-orbit gradient flow 
Eq. Q to the Lagrangian form is effected by taking 
Pghost = e~ 1 e^ + ir x e g e,g and taking the limit e — > 
to enforce the constraint 6 = d^x- Then the gradient flow 
becomes 

DO = _^6S 

Dt p(0 66' 

where /i(£) = 0(1) allows us to generalize the form as- 
sumed in Ref. [2] slightly, which we shall find necessary for 
reconciliation to be possible. 

5. Reconciling QFMin and Lagrangian Ghost 
surfaces 

We now seek to choose so that QFMin pseudo- 
orbits are also the Lagrangian ghost pseudo-orbits defined 
above. We require Q d^® = 0. That is, 

& = ®(6,0- (8) 

Then ¥ v = W, L v = L - v0, 8L y /d6 = dL/86, and mem- 
bers of our family of QFMin pseudo-orbits 

= m®<>) , (9) 
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where 0o is, as yet, an arbitrary label, satisfy the Euler- 
Lagrange equation 
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with© e (0,£) = d®(8,£)/d6. 

To reconcile QFMin and ghost orbits we require that 
the family of pseudo-orbits defined by Eqs Q and ( fT0] i is 
the same family as is generated by Eq. (JT]). Thus the la- 
bels ©o and r must be functionally dependent: t = t(0o), 
c/t = t'(0oV0o- Eliminating 6S/58 between Eq. ( 10 1 and 
Eq.Q and observing that ® g = (D®/D® )/(D8/D® ), 
where D8/D®o = d8(£\®o)/d®o, we find the reconcilia- 
tion condition 



D® 
D®a 



T'(0 () )v(0o) ID© 



D8 



(11) 



We now define ©o so that, for all ©0, 

r'(0 o )v(0„) = -1 , 
<9(<r|0 o + 2tt) = 6((\® )+2n, 



(12) 



choosing p(() so that Eq. (Hi satisfies the solvability con 



dition that the integral of both sides with respect to () over 
the interval [0, 2n\ must be 2n, giving 
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6. Perturbative construction of QFMin- ghost 
surfaces 

For example, consider 





L = — - e 2^ Vmji exp(z'm0 - in£) 



(14) 



with the reality condition V* = V_ m _„, and e the expan- 
sion parameter. (Defining the largest perturbation coeffi- 
cient, V2,i, to be 1, e = 10~ 3 for the case shown in Fig. [3]) 
As the unperturbed system is integrable, the expan- 
sions of v(@o), //(£), and ®(8, £) are of the form 

v = ev\ + e 2 V2 + . . . , 

= efii + e 2 ^ + ■ • ■ , (15) 
© = + 2(e©2, ) „ + ...)exp/(m0-«<r) . 

m.n 

and the (p, q) QFMin pseudo-orbits are expanded as 

0«1©o) = ^ + ©o + X(<» + e2 €i + ---) 

m,n 

x exp iUmtp^ - n) £ + m@o] ,(16) 



where t p>q = p/q. 

At 0(e) we find p\ = 0. Also, ©^„ is not used in the 
calculation of #„,„. Thus, to first order, uncorrected ghost 
and QFMin pseudo-orbits are identical. This explains why 
the low-order uncorrected ghost and QFMin orbits appear 
almost identical in Fig. [3] 



Explicit expressions for v\, ^2, (hnji> ^m,n an ^ ®mli ^ 
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given in Ref. 4, but, as there is a misprint in 9 mn , we give 
below a corrected version, 
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where N e m m , n n , = m'(m + m') \m + m' + m' '8( m +m')p,(n+n') q \ 
Also, 5 mpm = 1 - 6 mpm , where 8 mpm = 1 when mp = 
nq, is otherwise, and the prime on the sum over m' and 
«' indicates that the resonant terms, m'p = n'q, are to be 
deleted. 

7. ghost surface formulation 

In Sec. |4] we defined ghost surfaces in terms of an 
action-gradient flow in 9, but defined QFMin surfaces in 
terms of 0. Here we further unify the formulation by defin- 
ing the action-gradient flow also in terms of 0, 
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where T is the timelike label for ghost orbits [so dT = 
T'(&o)d@o] and the inverse function #(©, is defined by 
solving Eq. ([8]) for 9 as a function of 0. 



From Eq. ( 10 1 we also have, for corrected QFMin or- 
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Eliminating SS/S® between Eq. ( 18 i and Eq. ( 19 1 we find 
the new reconciliation condition 
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where, in the second equality, we have defined 0o so that 
r'(0o)v(0o) = -1, which satisfies the solvability condi- 
tion that C d®o D®/D®q - 2k without the need to intro- 
duce a factor like the required in Sec. [5] 

Furthermore, Eq. ( |20] > is satisfied by choosing the (- 
dependence of the pseudo-orbit labeling so that 



®((\®o) = O + <f 



(21) 



conjugating the 8 dynamics to straight-field-line dynamics 
(cf. Ref[6]). Thus the ghost surface approach appears to 
be much superior to 8 ghost surface approach, but needs to 
be tested in practice. 

As a first test we adapt the perturbative calculation in 
Sec.|6]to © ghost surfaces. First, expanding #(0, f ) = + 
Tm n 8m,n exp(zm0 - and using Eq. (21 



we recover 



Eq. \\6) . This ansatz is to be inserted in Eq. ( 19 1, written 
in the form 
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From Eq. ( 14 1, dLg /d£ = - Tm,n( mi ~ n) 2 8 m _„ exp[/m©o + 
i(mi- n)f], Lg = -e Tm.n ^ m ^m,n exp[/w#(£|©o) - inf], and 
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these into Eq. (22 1, and equating the LHS to zero at 0(e), 
we get 



VI (00) = -Yj^'np^VinJ"'^ , 



and 
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which are the same results as in Ref. 4, showing that 8 
and ghost orbits (and, by definition, reconciled QFMin 
orbits) are equivalent at this order. 
At 0(e 2 ) we find 
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and 
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where Nf l m> n n , = m'(m + m') [m + m' - m' 5 {m+m , )p , {n+ „, )q \ 
(the sign of the 6( m+m ') p ^„ + „')g term in this expression being 
the opposite of that for 9 ghost orbits). 

8. Conclusion 

We have reviewed the motivation and formulation of 
a recently published [4| unification of ghost and quadratic- 
flux-minimizing surfaces, and have presented a new for- 
mulation that appears more natural and more useful in that 
it could be used to form a straight pseudo-field-line coor- 
dinate system. 

More work remains to be done to develop a numeri- 
cal algorithm for constructing such ghost/QFMin surfaces 
at arbitrary order of perturbation, and to show they pre- 
serve the desirable properties of our previous formulation 
of ghost surfaces. 

The relation to heat transport in chaotic fields also 
needs to be investigated more deeply. 
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